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The one-particle reduced density matrix, p{x,y), of a many-electron quantum system 
can reveal a good deal about the presence or absence of spatial uniformity. For a general 
class of tight binding models, including the Hubbard model and the Falicov- Kimball model, 
P{x,y) turns out to have a particularly simple form in the case of the half- filled band. 
The result, which is applicable to a bipartite lattice, is that p{x,x) is always exactly 
equal to one for all a; on a finite lattice, even though the hopping matrix elements and 
interaction are nonuniform, random and uncorrelated. Furthermore P{x,y) = when 
X ^ y but X and y are in the same sublattice. This fact has been long appreciated 
in the chemistry literature, in the context of certain models describing tt electrons in 
conjugated carbon systems"*^. It was first observed for the Hiickel (free electron) model^ 
by Coulson and Rushbrooke^, who used it to justify the assumption that the effective 
potential should be the same at each carbon site in a self-consistent molecular orbital 
treatment. MacLachlan^'^ extended the result via a hole-particle symmetry argument to 
the Pariser-Parr-Pople (interacting electron) model^'^, of which the Hubbard model is 
formally a special case. For so-called alternant molecules, i.e., those in which the carbon 
atoms form a bipartite lattice, MacLachlan showed that to each N electron eigenstate 
corresponds a 2|A| — AT electron eigenstate with the complementary density. (Here |A| is 
the number of carbon sites, and there are 2 allowed spins at each site). These states have 
the same energy, up to a shift which is proportional to | A| — A". He used this pairing theorem 
to explain the identical spectra which had been observed for positive and negative ions of 
the same alternant molecule. Despite its usefulness, the result has remained unknown in 
the statistical mechanics and solid state physics literature, where the same models are used 
to investigate phase transitions and the existence of long range order. In this context, the 
persistence of uniformity in the face of randomness is striking. For instance, it hints at the 
stability of some periodic structures in solids, which is to say that the occurrence of periodic 
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structures in a system might be insensitive to some of the details of the Hamiltonian of the 
system. We therefore present the theorem in the context of statistical mechanical ensembles 
at positive and zero temperature, together with an extremely simple proof. Along the way, 
we extend the applicability of the theorem to include Hamiltonians with spin dependent 
hopping and spin-flip interactions. Thus consequences are extracted for some additional 
models of current physical interest, such as the Falicov- Kimball model. Finally, the infinite 
volume limit is discussed; if the infinite volume Gibbs state has non-constant density then 
it is not unique, and the theorem guarantees the existence of another Gibbs state having 
the complementary density. 

To establish some notation, we consider a finite graph (lattice) consisting of sites 
labeled by x, y, etc. and edges (or bonds) connecting certain pairs of sites. We assume 
that the graph is bipartite, i.e., the vertices can be divided into two disjoint subsets A 
and B such that there is no edge connecting x and y if x and y are both in A or both in 
B. The total number of sites in A,AorB is denoted by |A|, 1^41 or \B\. We assume that 
1^1 > \B\. We are given a hermitian |A| x |A| hopping matrix T with elements t^y = t*^. 
These elements are nonzero only if x and y are connected by an edge; thus the elements 
txy are zero whenever x and y are both in A or both in B. (In particular, txx = for 
all X.) Physically, T originates in overlap integrals and is real in the absence of magnetic 
fields that interact with the electron orbital motion. 

It is easy to see that the nonzero eigenvalues of T come in opposite pairs: for every 
eigenvalue A there is an eigenvalue —A. The two corresponding eigenvectors 0^ and 
are conjugates in the following sense: if (f)^ = (/^, g^) with being the ^-part of (f)^ and 
^f'^ being the S-part of (f)^ then 0"^^ = {f^, —g^). Alternatively, if V = ~ ^' 

the diagonal unitary matrix that multiplies by —1 on the B- sites, then VTV = —T. It 
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is possible that T has zero-modes; indeed if |^| > then T will have at least |^| — \B\ 
zero-modes and all these A = eigenvectors have the form (j) — (/, 0). 

Suppose now that we have a half-filled band, i.e., N = |A| electrons. By virtue of the 
two spin states for each electron we have that the ground state energy of Hq = ^ txyC^Cy 
is E = 2 Yli\<o ^- (Note: If |A| is odd there is at least one zero- mode, and so this formula is 
correct even in this case.) The ground state of Hq might be degenerate, however (because 
of zero- modes). We define the density matrix for spin a in the ground state to be 

A<0 A=0 

We see that Tr = J2x Pcri^i ^) = 1-^1/2, as it should, and that p^- agrees with the /3 — > oo 
limit of the positive temperature density matrix, defined in the grand canonical ensemble 

by 

P^^ix, y)^Yl ^\^y^\v) e"^V(l + e-^^). (2) 

A 

Note that we have used zero chemical potential which, by virtue of the A, —A symmetry, 
always yields |A| as the average particle number. If there are zero-modes the ground state, 
and the p^ in the ground state, will not be unique. Eq. (1) serves to fix p^. for our purposes. 
The Gibbs state is always unique for a finite volume. 

Another Gibbs state with which we shall be concerned is the canonical ensemble. The 
density matrix here will be denoted by PjSaix, y). Its definition is well known and we shall 
not write it explicitly for Hq, but we note that the /? — > oo limit of p^o- also equals p^-. 

Coulson and Rushbrooke's^ observation regarding (1) is the starting point of our 
further analysis, li x ^ A and y & A then, using the fact that (j)'^{x) = (j)~'^{x), we have 
that 

all A 
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since the form an orthonormal basis. A similar remark holds for x,y E B. Thus, 
Pa{x, x) — ^ for all a; e A and Pa{x, y) — for x,y G A or x,y G B. As we shall see from 
the following general theorem (by specializing to zero interaction) the same conclusion 
applies to Ppaix, y) and 7>i3a{x, y)- 

THE GENERALIZED HUBBARD MODEL 

The interacting system we shall be concerned with is the generalized Hubbard model 
defined by the Hamiltonian (with spin dependent hopping) 

H = ^^ txya<iaCya + Yl C^a;j,ar (2n^^ - 1) (2/1^^ - 1), (3) 

a x,y£A a,T x,y£A 

where Uxa = c\.a^xa ^^"^ with Uxyar real (but not necessarily of one sign). T^- = {txya} is 
hermitian and bipartite for each a =t or j. If we take — and U^yar — US^y then H 
is the usual Hubbard Hamiltonian (apart from a trivial additive constant) with interaction 
SUY^nx^Uxi- The noninteracting case, Hq, corresponds to Uxyar = 0. In general, the 
total spin angular momentum {SU (2) symmetry) will be conserved if we require txya and 
Uxyar to be independent of the spin labels a and r; for our purposes we do not require 
this SU (2) invariance. 

The positive temperature, grand canonical density matrix p^^- is defined to be 

p^a{x,y) = Z-'Tr[cl,Cy,e-^% (4) 

where Z = Tr [e~^^] . Formula (4) reduces to (2) for the noninteracting case. The trace 
is over the full Fock space containing all particle numbers ranging from to 2|A|. Again, 
the zero chemical potential in (4) insures that Tr p^^r = |A|/2. 

The canonical density matrix Pi3a{x,y) for this model is also given by (4), but where 
the trace is only over the A'"-particle sector (note that both H and c^^Cy^ leave this sector 
invariant). The half-filled band is defined by iV = |A|. Since {Tr ABCy = Tr C^SUt = 
Tr B^A^C''' we see that Pi3a and p^a are hermitian matrices for each a. 
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THEOREM: (Uniform Density in the Generalized Hubbard Model). 

The canonical and the grand canonical density matrices satisfy: 

Pl3a{x,x) = Pi3aix,x) = 1/2 for all x G A (5) 
Pl3a{x,y) = P(3a{x,y) = H X, y ^ A OT X, y ^ B . (6) 

Proof: The proof for PjSa- will be the same as that for PjSa- so we shall only give the proof 
for the latter. 

First, we consider real T. MacLachlan's version^ of the hole-particle unitary transfor- 
mation, 

Cxa ^ cla for X eA, c^^ ^ -c\^ for x G S, (7) 

evidently leaves the Hamiltonian H and the relevant Hilbert spaces invariant. If this 
unitary transformation is denoted by W we have that W'^ = 1, = W'^ ^ WHW = H and 
hence ZPf3^{x,y) = Tr[Wcl^Cy^WWe-f^^W] = Tr[{Wcl^W){Wcy^W)e-'^^]. If x,y e A 
we can use (7) and the fermion commutation rule to conclude that Pi3a{x,y) = Sxy — 
Pf3a{y, x). The same is true if x,y e B. If T is real, ppa is evidently real; since P/Sa is also 
hermitian the theorem is proved in the real case. 

The complex case is a bit subtle. The Hamiltonian H is no longer invariant under the 
hole-particle transformation W, but it is invariant under the antiunitary transformation 
Y = JW, in which J is complex conjugation. More precisely, any vector ^ in our Hilbert 
space can be written as a linear combination, with complex coefficients, of the basis vectors 
consisting of monomials in the cj,^'s applied to the vacuum. The antiunitary map J acts 
on by replacing each coefficient by its complex conjugate. We note that JW = WJ and 
therefore = 1. It is also easy to see that YHY = H and that YcxaY = Wcx<jW, which 
is given by (7). 
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Now suppose that K is an arbitrary linear operator, and consider L = JKJ. Although 
J is nonlinear, it is easy to check that L is linear. In fact the matrix elements of L in the 
above mentioned basis are simply the complex conjugates of the corresponding elements of 
K. Therefore, even though = 1, it is not generally true that TrL = Tr JKJ = Tr KJ^ = 
TtK. What is true is that TrL = (TrK)*. 

In our case we have, for x,y E A or x,y E B, Zpp^{x, y)* = Tr[JWcl.^Cy^e~^^W J] 
= Tr[{YcUY){YCy,Y){Ye-^^Y)] = IY[c,,4,e-^^] = Z {5,y - Pp,{y, x)). The hermitic- 
ity of p^o- now implies the theorem. QED. 

It is worth noting that only the invariance of H under Y and the bipartite structure 
of the lattice have been used. Thus the theorem (with some obvious modifications) applies 
to the case where the kinetic energy has spin-flip terms and is given by 

x,y,a,T 

Hermiticity requires t^yar = tlxra'i be a hopping model on a bipartite lattice imposes 
another condition on T which will be made clear in a moment. To apply the previous 
theorem, we observe that this model is equivalent to a system of spinless fermions living on 
a lattice twice as large as A, with new co-ordinates (x, a), where x E A and a e (j, J,}. The 
new lattice should be bipartite, which is to say that it can be divided into two sublattices 
such that txyar — when both {x,a) and (y, r) are in the same sublattice. Pi3a{x,y) 
is replaced by Pi3{x,a;y,T) = Pi3{y,T;x,a)*, and the extended theorem states that it 
equals ^Sj^yS^r when {x,a) and (y,T) are both on the same sublattice. In practice, the 
only subtlety is the identification of the sublattices. For example, let A be decomposed 
into A and B. Then the theorem applies if txyar = whenever x,y e A or x,y e B. 
In this case one of the sublattices would consist of those sites (a, a) for which a E A 
and cr E {T; i}- Alternatively, we could have taken the condition to be txyaa = when 
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x,y & A OT x,y & B, and txy^i = whenever one of x,y is in A and the other in B. 
This model also satisfies the hypothesis, but in this case one of the sublattices consists of 
{ (a,^) \ a E A} [J { {b, I) \ b E B }. The latter scenario would be realized if we had only 
on-site spin-fiip terms: tx^tJ. ~ tx^xy, a physically appealing possibility. In the interest of 
simplicity of notation and exposition we have relegated this generalization — relevant to 
some physical models — to the remark here. 

The theorem also extends to include spin-spin interactions of a more general nature 
than those manifested in (3). This is particularly welcome, since many physical models 
— the t-J model for example — require SU(2) invariant spin couplings. Let x E A and 
consider the operators = {uxi — nxi)/2 and = c]^-|-c^|- They generate the usual 
SU(2) algebra at site x: 

= ( {S+ + S-)/2, {S+ - S-)/2t, Si) 

where S~ = (S^)"^ and the total spin at x is Sx ■ Sx- In terms of these operators, the 
spin-spin interaction realized in (3) is of the form S^Sy where x,y E A; it occurs with a 
coupling constant proportional to {Uxy-f-i + Uxyii) — {Uxyi i + Uxyi-i) ■ Since S^. ^ —S^. under 
the hole-particle inversion Y, the invariance of this interaction is manifest, but it clearly 
destroys the SU(2) symmetry of the Hamiltonian. However, the antiunitary transformation 
Y does not single out the 3-direction; indeed it follows that Sx ^ —Sx under Y. As a 
result, the conclusions of the uniform density theorem, (5) and (6), are valid for any system 
described by a Hamiltonian 

3 

H' = J211 txyAaCya + E ^-^(^^ " " 1) + E E '^S^^'^^- 

cr x,y€A x,y€A x^y i,j=l 

Here = Uxi + riyi and {S^, S^, S^.) — Sx, txya is as in (3) and the Uxy and J*-^ = 
1,2,3) are completely arbitrary real constants. The case J*-^ = Jxy^ij corresponds to the 
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SU(2) invariant interaction S^- Sy. We have omitted the spin-charge couphngs (ux — l)Sy, 
which would appear in the most general form of H' , in order to maintain its physical 
simplicity. 

We conclude by discussing some consequences and extensions of the Theorem in certain 
limiting cases. 

I. No Correlations within a Sublattice for the Free Electron Model: 

For the free electron model, the interaction terms Uxyar are zero, leaving a purely 
kinetic Hamiltonian of the form (3) — or (8), although we continue to suppress this 
generalization for notational convenience. This Hamiltonian is quadratic in the creation 
and annihilation operators, i.e., it is a one-body operator. Wick's theorem® therefore 
applies to the grand canonical ensemble^, for which it states that all operator product 
expectations are expressible in terms of pair expectations. We use this observation together 
with the vanishing of Ppai^i v) for a;, y in the same sublattice to show that there can be no 
correlations between the electron densities on sites in this sublattice. This striking result 
does not seem to appear previously in the literature, although it generalizes a comment 
made on p. 86 of Salem's book^ about two particle correlations in the ground state. 

Electron correlations are manifested in the n-particle reduced density matrices p^'^^ 
which, for the grand canonical ensemble at inverse temperature are defined by 

P^pa^^'^i ■ ■ ■ 1 ^n'l yi, ■ ■ - Vn) = Z ^Tr [cj-j^cj-^o- • • • cl^a'^y^a- ■ ■ ■ ^yicr'^yiCT^ • 

The trace is over the full Fock space containing all particle numbers to 2|A|, and the 
normalization corresponds to Tr p^^} = Tr[N^{N^ -!)■■■ {N^ -n + l)e-^^]/Z, where 
counts the number of spin a particles. 
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Diagonalizing the hopping matrix through a change of basis, the Hamiltonian becomes 
H — ^x^^\^x- Here A runs over 2|A| values. b\ and are quasi-particle creation and 
annihilation operators, and are related to cl^ and c^,^ through a unitary matrix 0: 



A A 

Consequently, and 6^ obey the fermion anticommutation relations 



(9) 



{bib,} = 5,,, {blbi} = {bx,b^}=0. 

Of course, p^^J can be expressed in terms of 6| and b^ via (9), and it is in this form that 
it will be most easily evaluated. 

Letting nx = b^b^, the anticommutation relations imply 



Tr 



bl ■ ■ ■ blb,^ ■ ■ ■ b,,e-^''] = det [<Ja.mJ ■ [nx, ■ ■ ■ n^^e"^^] (10) 



since the left hand side vanishes unless the Aj are distinct and Hi = XTr{i) for some per- 
mutation TT e (Sn on n labels. The antisymmetry is manifested in the determinant. For 
the grand canonical ensemble, the trace is over the full Fock space, which is an (antisym- 
metrized) tensor product of the one particle Hilbert space. Taking the trace in the bx 
basis, H is diagonal so 

n 

Z-^Tr[nx,---nx^e-^'']=ll{l + e^'T'- 

1=1 

Because this is a product, it can be absorbed into the determinant in (10). The right 
hand side of (10) then becomes Zdet [Sxi^jil + e^^^)~^]^ As a particular instance, 

5A^(l + e^^)-i,thus, 



bib^e-^^ 



b{ ■■■b\ b^ ■■■b^e-f^^ 



= det 



Z-^Tr 



-pH 
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Multiplying both sides by (f)^^ ixi,a)* ■ ■ ■ 0^^" (xn, a)*(j)^^ (yi, u) • • • 0^" (yn, cr) and summing 
over all A's and //'s brings us back to the position basis through (9). Since the determinant 
is termwise multilinear in the h\, and h^^. , we obtain Wick's Theorem: 

P^p}{{xi^yi}i) = det 

Now take all xi, . . . , Xn, yi, ■ ■ ■ , i/n in the same sublattice, either A or B. By the Uni- 
form Density Theorem, 

Explicitly, P^^J{xi, . . . ,Xn',yi, ■ ■ -yn) = ±2~" when the Xi are distinct and Xi = for 

(n) 

some fixed permutation tt e Sn- The sign of Pj^J is selected by the parity of tt. Otherwise, 

= 0. 

This is to say that there can be no spatial correlations in the electron density between 
sites of the same sublattice, save only that the probability of finding more than one electron 
in the same site and spin state vanishes. As an example, if x,y E A (or B) the cr-clectron 
pair density is P^pl{x,y;x,y) = (1 — Sxy)/4:, the product of the one-particle densities at 
X ^ y. Thus the Pauli pressure is not felt between the electrons at these sites. This is 
surprising since one generally expects, and indeed will find, correlations between arbitrary 
sites X E A and y E B. 

The extension to Hamiltonians of the form (8) is immediate. The n-particle reduced 
density matrices now depend on spin as well as spatial co-ordinates, and are defined in the 
obvious way. Taking all the {xi, (Ji) and (i/i, r^) to be in the same sublattice, the conclusion 
is that 

P^p^ (xi, CTi, . . . , Xn, CTn; 2/1, Ti, . . . , 2/n, Tn) = 2~"det [^a;ij//aiT,- ] ^^-^^ • 



11 



EHLMLRJM-22/Peb/93 

Unfortunately, these results extend neither to the canonical ensemble nor to the inter- 
acting Hamiltonian. For the canonical ensemble, correlations result from the fixed number 
of available particles. To see this, consider a six site lattice in which four of the sites 
are totally isolated, and the remaining two are joined by a a non-zero hopping matrix 
element t. The hopping matrix eigenvalues are with multiplicity four, and ±|t|. If this 
lattice is populated with spinless fermions, then each of the 16 ground states will involve 
one particle hopping between the connected sites, and zero to four particles distributed 
over the isolated sites. However, only 6 of these states, those having a total of three 
particles, will be represented in the canonical Gibbs state. Select any two of the isolated 
sites, and consider the probability of finding them both occupied. In the grand canonical 
ensemble, we have seen that this must be 1/4 at any temperature. However, in the zero 



temperature canonical ensemble, the likelihood of this event is 




= 1/6. Of course. 



this is a finite size effect. On the other hand, it is easy to create transparent counterexam- 
ples for the interacting case by taking T = 0. In particular, the Ising model Hamiltonian 
H = JijSiSj can be embedded in the interaction term of (3). The identifications 

are L = A x {|, J,}, Sj = Sn^jcr — 1 where i — {x,a), and J^- = Uxyar so that an occu- 
pied spin-site state corresponds to an up Ising spin. This model is well known to have 
correlations between all sites for a connected lattice. 

Actually, we have just discussed two trivial limits for the standard Hubbard model 
in which Uxyar = US^y- The U ^ limit is the free electron model, while C/ — > oo (or 
equivalently T — > 0) corresponds to a (spatially) disconnected Ising model. In both of 
these limits all spatial correlations which are detectable on a sublattice must vanish — at 
least for finite lattices and (3 < oo, where the observables are continuous functions of U. In 
particular, the antiferromagnetic correlations (which can be computed to be present on a 
small lattice, and are conjectured to be long range for £) > 3 dimensional square lattices) 
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are clearly seen to be a result of interference between the kinetic and potential terms in 
the Hamiltonian. 

II. The Falicov-Kimball model^^: In this model there is one species of mobile, 
spinless electrons and one species of arbitrarily fixed particles. It is the same as the 
Hubbard model with the choice txyi = and txy^ — txy The Uniform Density Theorem 
applies to this model as a special case. Note that it says that both the mobile and immobile 
particles have density 1/2. 

III. Ground States: If we define the ground state Pa as the limit /3 — >^ oo of the 
canonical Pp^-, then (5) and (6) apply there, too. (Note: For the interacting Hamiltonian, 
it is not clear that Pp = P/3 in the /3 — > 00 limit. Certainly the canonical ensemble and the 
grand canonical ensemble do not coincide in this limit, even for the free electron model. 

In Remark I. above, we gave a contrived example involving a six site lattice for which 

(2) ~ (2) 

Pp ' at /3 ^ 00, in the process of showing that certain correlations did not vanish in 
the canonical Gibbs state.) Of course, the ground state is not generally unique, and there 
are other possibilities for p^-) in which case (5) and (6) apply to states that are invariant 
under Y = WJ. Not every ground state is Y invariant. In the case of the usual Hubbard 
model on a connected lattice with U > and real Y it is known^^ that the ground state 
has spin angular momentum S = {\A\ — \B\)/2 and it is unique apart from the (25' + l)-fold 
degeneracy associated with — |(-^T ~-^J.) ^ —5' + 1, . . . , +S}. Ground states that 
are Y invariant are mixtures of states with and —S^, i.e., ^\S^){S^\ + || —S^){—S^\ in 
Dirac's notation. The reason we can be sure that Y\S^) = \ —S^) is this: Y\S^) is — in 
any event — a state with |(iV-|- — N^) — —S^. It is also a ground state. By uniqueness, 
this state must be | —S^). In the general model (3), we cannot be sure that spin flip = 
hole-particle transformation. 
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IV. Infinite Volume Gibbs States: These states can, of course, have different 
properties from finite volume states. One way to define them is as limits of finite volume 
states with specified boundary conditions that need not respect Y symmetry. One example 
concerns the Falicov-Kimball model with t/ > on a hypercubic lattice in D dimensions: 
For D > 2 it has long-range order in the ground state and at low temperatures, in which 
the up spins preferentially occupy the A-sites and the down spins the 5-sites (or vice 
versa) Similarly, the usual Hubbard model is expected to show the same behavior 
when D > 3 ii U > 0, at least if t/ > is large enough. This has not been proved, 
however. These examples violate (5), but they do suggest that the charge density satisfies 
Pl3^{x,x) + Pi3i{x,x) = 1 in these models. If we now introduce nearest neighbor repulsion 
(which is allowed in our general model) then even this constancy of the charge density 
might be violated, however. 

What is significant about our finite volume theorem (5), (6) — and which does remain 
true in the infinite volume limit — is that for every state with non-constant (spin or 
charge) density there is another Gibbs state corresponding to the same temperature which 
has the complementary density. In other words, one cannot invent a non-translationally 
invariant Hamiltonian (either by altering the hopping matrix or the pair potentials) with 
the property that it forces the density to increase in some specified regions and to decrease 
in others — even though one might have thought a-priori that the density can be controlled 
by the hopping or potential energy. Any attempt to cause a non-constant density will 
always result in the certainty that exactly the reverse of the desired non-constancy will 
occur with the hole-particle reversed boundary condition. 
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